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Abstract. A measurement of the virtual-photon asymmetry A2{x,Q^) and of the spin-structure function 
g2{x,Q^) of the proton are presented for the kinematic range 0.004 < j; < 0.9 and 0.18 GeV^ < < 20 
GeV^ The data were collected by the HERMES experiment at the HERA storage ring at DESY while 
studying inclusive deep-inelastic scattering of 27.6 GeV longitudinally polarized leptons off a transversely 
polarized hydrogen gas target. The results are consistent with previous experimental data from CERN 
and SLAC. For the a;-range covered, the measured integral of 52(2;) converges to the null result of the 
Burkhardt-Cottingham sum rule. The moment of the twist-3 contribution to 52(2;) is found to be 
compatible with zero. 
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The description of inclusive deep-inelastic scattering 
of longitudinally polarized charged leptons off polarized 
nucleons requires two nucleon spin-structure functions, 
gi{x,Q^) and g2{x,Q^), in addition to the well-known 
structure functions Fi{x,Q^) and ^2(2;, Q^) [I]. Here, 
—Q^ is the squared four- momentum of the exchanged vir- 
tual photon with laboratory energy v, x = Q^/{2Mv) 
is the Bjorken scaling variable, and M is the nucleon 
mass. In the quark-parton model (QPM), the spin struc- 
ture function gi^x^Q"^) can be interpreted as a charge- 
weighted sum of the quark helicity distributions Aq(x, Q^) 
describing a longitudinally polarized nucleon. 



9i{x, 



(1) 



The spin structure function g2{x,Q^) does not have 
such a probabilistic interpretation in the QPM. Its prop- 
erties can be interpreted in the framework of the operator 
product expansion (OPE) analysis [2ll3|4], which shows 
that g2{x,Q^) is related to the matrix elements of both 
twist-2 and twist-3 operators. Neglecting quark mass ef- 
fects, g2{x, Q^) can be written as a sum of two terms 

g2{x,Q') = g^'^{x,Q') + -g2{x,Q') . (2) 

Here, g^^ {x, Q^) is the twist-2 part derived by Wandzura 
and Wilczek [S] 
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'^'^{x,Q') = -g,{x,Q^) + 
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(3) 



The second term in Eq. g2{x, Q^), is the twist-3 part 
of (72(2;, Q^)- It arises from quark-gluon correlations in the 
nucleon and is the most interesting part of the function. 
The moment oi g2{x,Q^), 



rf2(Q') = 3 / x^Mx, 
Jo 



Q^)dx, 



(4) 



can be calculated on the lattice (see, e.g., [HIIZ], where ^2 
is defined with an additional factor of two with respect 
to (HI). The moment 6.2 has also been linked to the trans- 
verse force acting on the quark that absorbed the virtual 
photon in a transversely polarized nucleon, and thus to 
the Sivers effect Ii[9l[10|. 

The Burkhardt-Cottingham sum rule [TT] for g2 at 
large Q^, 



/ g2{x,Q'')dx = Q, 
Jo 



(5) 



does not follow from the OPE. Its validity relies on an 
assumed Rcgge behaviour of 92 at low x. In the absence of 
higher twist contributions to the function g2, i.e., g2{x) = 
0, the sum rule would automatically be fulfilled. Hence a 
violation of the sum rule would indicate the presence of 
higher- twist contributions. 

The spin structure functions gi{x,Q^) and 32(2;, Q^) 
can be related to the virtual photon-absorption asymme- 
tries Ai{x,Q'^) and A2{x,Q'^) [J 



A. = 
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Here, cr^2 ^-'^'^ '^T/2 ^''-'^ transverse virtual-photon ab- 
sorption cross sections for total photon plus nucleon an- 
gular momentum projection on the photon direction of 
1/2 and 3/2, respectively. The cross-section a^^ arises 
from the interference between the transverse and longitu- 
dinal photon- nucleon amplitudes, with 7 = 2Mxj \/ Q^. 
All of the cr's are differential cross sections depending on 
X and Q^, but this dependence was omitted for brevity. 

The measurement of the structure function 32 requires 
a longitudinally polarized beam and a transversely po- 
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larized target. In this case, the inclusive differential cross 
section can be represented as a sum of two terms, the 
polarization-averaged part, (T[/[/, and the polarization- 
related part, ctlt- Here, the subscript uu indicates that 
both the beam and the target are unpolarized, while the 
subscript lt indicates a longitudinally polarized beam 
and a transversely polarized target. The polarization- 
related part of the cross section at Born level, i.e., in 
the one-photon approximation, is given by [3] 



with 



dxdydcj) 



hi cos ( 



i-y 



(|.9i(a:,Q') + <72(a;,Q2)) 



(8) 



Here, hi = +1 (—1) for a lepton beam with positive 
(negative) helicity, a is the fine-structure constant, and 
y ~ v/E, where E is the incident lepton energy. The an- 
gle (f> is the azimuthal angle about the beam direction be- 
tween the lepton scattering plane and the "upwards" tar- 
get spin direction. The polarization-related cross section 
aLT is significantly smaller than the polarization-avera- 
ged part auu and therefore its measurement requires high 
statistical precision. Up to now, the function 52 and the 
asymmetry A2 have been extracted [T^[T51[n] to less ac- 
curacy than gi and Ai. 

A measurement of the inclusive cross sections ([5]) at 
angles cj) and 4> + Tr allows one to construct the asymmetry 
Alt, 



A^,T(x,g^0) = h, 



a{x, Q'^,4>) - a{x, Q'^,(j) + tt) 
(t{x, Q'^,(j)) + a{x, -f tt) 



^h, 



(Tlt{x,Q^ 

(7uu{x,Q'^ 



'At{x,Q^) cose/}, (9) 



which defines the asymmetry amplitude At{x,Q'^). This 
amplitude contains all information on the function 52 and 
the asymmetry A2. Their extraction requires the knowl- 
edge of (Tuu{x, Q^), which can be expressed by the struc- 
ture functions i^i_2(x, Q^) or, equivalently, parameteriza- 
tions of the function F2{x, Q^) and the ratio of longitudi- 
nal to transverse virtual-photon absorption cross sections 
R = R{x,Q^). The extraction of the structure function 
g2{x,Q^) from the asymmetry amplitude At is analo- 
gous to the extraction of gi{x, Q"^) from the longitudinal 
asymmetry as described in |15| . The function (72 can be 
extracted from the measured asymmetry amplitude At 
and parameterizations of previous measurements of auu 
and gi, using ^ and Also Fi can be computed from 
parameterizations of F2 and R. This leads with ([7]) to the 
following relations 



Pi f^T , gi 



A2 = 



1 



1 + 7^ 
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y{2-y){l + j^y/2) 



y2(l+^2) + 2(l-2;-7V/4)(l + i?) 



(12) 
(13) 

(14) 



However, it is not obvious from these relations that the 
extraction of g2 is independent of correlated variations in 
values of Fi, F2 and R that conserve the directly mea- 
sured values of auu ■ 

This paper reports a new measurement of the func- 
tion 52 and the asymmetry A2. The data were collected 
during the years 2003 - 2005 with the HERMES spec- 
trometer [TH] using a longitudinally polarized positron or 
electron beam of energy 27.6 GeV scattered off a trans- 
versely polarized target [17j of pure hydrogen gas inter- 
nal to the HERA lepton storage ring at DESY. The us- 
age of a pure target avoids the complications of nuclear 
corrections present in previous measurements. The open- 
ended target cell was fed by an atomic-beam source [18] 
based on Stern-Gerlach separation combined with radio- 
frequency transitions between hydrogen hyperfine states. 
The nuclear polarization of the atoms was flipped at 1-3 
minute time intervals, while both the polarization mag- 
nitude and the atomic fraction inside the target cell were 
continuously measured [19j . The average magnitude of 
the proton polarization was 0.78 ±0.04. The lepton beam 
(positrons during 2003 - 2004 and electrons in 2005) was 
self-polarized in the transverse direction due to the asym- 
metry in the emission of synchrotron radiation |20j in the 
arcs of the HERA storage ring. Longitudinal orientation 
of the beam polarization was obtained by using a pair 
of spin rotators [H] located before and after the inter- 
action region of the HERMES spectrometer. The sign 
of the beam polarization was reversed every few months. 
The beam polarization was measured by two independent 
HERA polarimeters p2l[23l[24] . The average magnitude of 
the beam polarization was found to be 0.34 ± 0.01. The 
scattered leptons were detected by the HERMES spec- 
trometer within an angular acceptance of ±170 mrad 
horizontally and ±(40 — 140) mrad vertically. The lep- 
tons were identified using the information from an elec- 
tromagnetic calorimeter, a transition-radiation detector, 
a preshower scintillating counter and a dual-radiator ring- 
imaging Cerenkov detector. The identification efficiency 
for leptons with momentum larger than 2.5 GeV exceeds 
98%, while the hadron contamination in the lepton sam- 
ple is found to be less than 1%. The luminosity monitor 
pS] measured e+e~ (e~e~) pairs from Bhabha (M0ller) 
scattering of beam positrons (electrons) off the target gas 
electrons, and 77 pairs from e+e~ annihilation in two 
NaBi(W04)2 electromagnetic calorimeters, which were 
mounted symmetrically on either side of the beam line. 
Tracking corrections were applied for the deflections of 
the scattered particles caused by the vertical 0.3 T target 
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holding field, with little effect on the extracted asymme- 
tries. 

Most of the details of the analysis follow the inclusive 
analysis described in |15j . The kinematic constraints im- 
posed on the events were: 0.18 GeV^ < Q'^ < 20 GeV^, 
invariant mass of the virtual photon-nucleon system W > 
1.8 GeV, 0.004 < x < 0.9, and 0.10 < y < 0.91. After 
applying data quality criteria, 10.2 x 10^ events were 
available for the asymmetry analysis. The kinematic re- 
gion covered by the experiment in {x, (3^)-space was di- 
vided into nine bins in x. Each of the seven x-bins in the 
region x > 0.023 was subdivided into three logarithmi- 
cally equidistant bins in Q^. The range in 0-space (27r) 
was divided into 10 bins. Two of the 0-bins cover the 
shielding steel-plate region of the spectrometer and thus 
cannot be used for the analysis. The data were corrected 
for the e'^e~ charge-symmetric background |15| . which 
amounted in total to about 1.8% of the events, reaching 
the largest contribution of about 14% at small values of 

X. 

The measurement of the asymmetry ALT{x,Q^j(f>) 
given by ^ can be performed by either reversing the 
transverse target polarization and comparing the num- 
ber of events in the same part of the detector, or by com- 
paring the number of events in the upper and lower part 
of the detector for the same upward or downward target 
polarization direction. The first method provides a bet- 
ter cancellation of acceptance effects and was chosen to 
obtain the asymmetry 

ALrix.Q"^ ,cj),hi) = 

' mii^{x,Q'^,c^)Cp'^ + Nhi^{x,Q-^,(t>)Cp'^ ■ 

Here, N^^^^^"^ is the number of scattered leptons in one 
bin of the 3-dimensional space {x,Q'^,(j)) for the case of 
the incident lepton with helicity hi when the direction 
of the proton spin points up (down), ^'^'l^t^) and c'p^''^^ 
are the corresponding integrated luminosities and the in- 
tegrated luminosities weighted with the absolute value of 
the beam and target polarization product, respectively 

^ Jdt (t)r(t) , (16) 

^lum) ^ J dt L'"^'^'i-\t) I PB{t)PT{t) I T{t). (17) 

Here, L{t) is the luminosity, T{t) is the trigger live-time 
factor, and Pb and Pt are the beam and target polariza- 
tions, respectively. The asymmetries evaluated according 
to were found to be consistent for the two beam 
helicity states. Therefore they were combined in the fur- 
ther analysis. Finally, the asymmetry given by ([T5|) was 
unfolded for radiative and instrumental smearing effects 
to obtain the asymmetry corresponding to single-photon 
exchange in the scattering process. Radiative corrections 
were calculated using a Monte-Carlo generator [26]. The 
unfolding procedure is analogous to that used previously 
in other HERMES analyses [T51[?7l[^ . It inflates the size 



of the statistical uncertainties especially in the lowest Q^- 
bins at a given value of x. The magnitude of inflation 
reaches almost a factor of two at low values of x. The 
subdivision of x-bins in the range x > 0.023 into three 
bins in decreases the error inflation by about a fac- 
tor of 1.5 because at larger Q"^ the amount of smearing 
between x-bins is smaller and the prefactors of At in 
(jlOp and pT|) are larger in magnitude. After the unfold- 
ing procedure the central values of g2 and A2 changed 
less than the initial statistical uncertainties. As a con- 
sequence of the unfolding procedure, the resulting data 
points are no longer correlated systematically through 
radiative and instrumental smearing effects, but are only 
statistically correlated [TS]. The procedure generates a 
statistical covariance matrix for the data points. 

In every {x, Q^)-bin the amplitude At{x, Q^) was ob- 
tained by fitting the unfolded asymmetries with the func- 
tion /(0) = —At{x,Q^) cos<j>. Finally, the asymmetry 
A2{x, Q^) and the function g2{x, Q^) were evaluated from 
the amplitude At and the previously measured function 
gi , for which a world-data parameterization |29| was em- 
ployed, using ([TU)) and pT|) . The structure function 

Fi(x, Q2) = F2{x, Q2)(i + ^^)/[2x{l + R{x, Q^))] (18) 

was calculated using a parameterization of the structure 
function F2{x,Q'^) [3D] and the ratio R{x,Q'^) [SI]- AU 
kinematic factors in ([TUj) and ([TT|l . and the functions Fi 
and gi / Fi were calculated at the average values of x and 

in each (a;, Q^)-hm after unfolding. 

The uncertainties in the measurements of the beam 
and target polarizations produce in total a 10% scale un- 
certainty on the value of At- Other sources of system- 
atic uncertainties such as acceptance effects, small beam 
and spectrometer misalignments, uncertainties in the tar- 
get polarization direction, correction for track deflection 
in the vertical target holding field, the unfolding proce- 
dure and a possible correlation between prefactors of At 
and At itself in (fTPj) and (|TT]) were evaluated by Monte- 
Carlo studies. Uncertainties stemming from parameteri- 
zations of gi{x,Q^), F2{x,Q^), and R{x,Q'^) were esti- 
mated also. In the error propagation to g2, the uncer- 
tainty in R{x, Q^) was not included in addition to that of 
F2(x, Q^), since they are strongly correlated as explained 
in jl5j . The total systematic uncertainty was evaluated 
as the quadratic sum of all the considered sources. Its 
magnitude is less than the magnitude of the statistical 
uncertainty. 

Figure[T]shows the values of xg2 as a function of for 
the bins with x > 0.1, which have sufficient coverage in 
along with results from the E143 [13] and E155 [14] 
experiments at SLAC. The entire set of measured data 
and average values of x and are presented in Table [1] 
Within the accuracy of the data, they are in agreement 
with the other experiments. Also shown is the Wandzura- 
Wilczek term 32^^, which was evaluated according to 
A world data parameterization oi gi{x, Q^) [25] was used 
for the calculation. 

In order to study the x dependence, A2{x,Q^) and 
g2{x, Q^) in bins covering the same x range but with dif- 
ferent values were evolved to their mean value of 
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Fig. 1. The spin-structure function xg2{x, Q^) of tiie proton 
as a function of for selected values of x. Data from the 
experiments E155 [14] and E143 [13] are presented also. The 
average values of x for these two experiments are slightly dif- 
ferent from the HERMES values of (x) indicated in the panels. 
The error bars represent the quadratic sum of the statistical 
and systematic uncertainties. The solid curve is the result of 
the Wandzura-Wilczek relation ((3} 



and then averaged. The evolution of A2{x,Q'^) was car- 
ried out assuming that the product \fQ^A2 does not 
depend on Q^, which follows from ([7]), since gi/Fi is 
known to vary only weakly over . The structure func- 
tion g2{x,Q^) was evolved assuming that its depen- 
dence is analogous to that for the Wandzura-Wilczek part 

of 52- 

The averaged results for xg2 and A2 and the statis- 
tical and systematic uncertainties are listed in Table [2l 
where the average values of x and are also given. The 
quoted statistical uncertainties correspond to the diago- 
nal elements of the covariance matrix obtained from the 
unfolding algorithm. The correlation matrix for xg2 in 
nine .x-bins is presented in Table [3P1 

The results for the virtual-photon asymmetry A2 and 
the spin-structure function xg2 as a function of x arc pre- 
sented in Fig. [2] together with data from the experiments 
E155 [g, E143 [HI, and SMC [Hj. The HERMES data 
are shown for two regions of Q^, (Q^) > 1 GeV'^ (closed 
symbols) and {Q^) < 1 GcV^ (open symbols). The exper- 



It is also available in 23 bins for the data in Table [T] 
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Fig. 2. Upper panel: The virtual-photon asymmetry A2 
of the proton as a function of x. Bottom panel: The spin- 
structure function xg2 of the proton as a function of x. HER- 
MES data are shown together with data from the E155 |14) . 
E143 [13], and SMC [TJ] experiments. The total error bars 
for the HERMES, E155, and E143 experiments represent the 
quadratic sum of the statistical and systematic uncertainties. 
The statistical uncertainties are indicated by the inner error 
bars. The error bars for the SMC experiment represent the sta- 
tistical uncertainties only. The solid curve corresponds to the 
Wandzura-Wilczek relation ^ evaluated at the average 
values of HERMES at each value of x. For the HERMES data, 
the closed (open) symbols represent data with (Q'^) > 1 GeV'^ 
({Q2) < 1 GeV^) 
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iments have only slightly different values of average 
for a particular value of x. The results are within their un- 
certainties in good agreement with each other. The solid 
curves represent values of A2 and xg2 evaluated with the 
Wandzura-Wilczek relation ([3]) using the pa- 
rameterization The values were calculated at the 
average of HERMES at each value of x. Within the 
uncertainties the data satisfy the positivity bound |32j 
for the asymmetry A2, IA2I < \/ Ri\ + Ai)/2 ~ 0.4, for 
all values of x in the kinematic range of the HERMES 
experiment. 

The Burkhardt-Cottingham integral ([5]) was evalu- 
ated in the measured region of 0.023 < a; < 0.9 at = 
5 GeV^, resuhing in lH^^ g^ix, Q^) dx = 0.006 ± 0.024 ± 
0.017. This result is to be compared with the combined 
result from experiments E143 and E155 [T3| in the region 
0.02 < X < 0.8: /°g2 52(3;, Q2) dx = -0.042 ± 0.008. 

Using the results measured by HERMES for the func- 
tion g2, the twist-3 matrix element d2 given by (|4]) was 
evaluated. For the unmeasured region 0.9 < a; < 1, the 
ansatz 32(2;) oc (1 — x)^ was assumed. The uncertainty 
in the extrapolated contribution was taken to be equal 
to the contribution itself. The contribution from the re- 
gion X < 0.023 was neglected because of the x^ suppres- 
sion factor. The result is ^2 = 0.0148 ± 0.0096(sfai) ± 
0.004:8{syst). This is to be compared with the combined 
result from experiments El 43 and E155 (i2 = 0.0032± 
0.0017. 

In conclusion, HERMES measured the spin-structure 
function g2 and the virtual-photon asymmetry A2 of the 
proton in the kinematic range 0.004 < a; < 0.9 and 0.18 
GeV^ < < 20 GeV^. For the covered x-range the mea- 
sured integral of 172(2;) converges to the null result of the 
Burkhardt-Cottingham sum rule. The moment of the 
twist-3 contribution to g2{x) is found to be compatible 
with zero, in agreement with expectations on its small- 
ncss from lattice calculations. The results on A2 and g2 
are overall in good agreement with measurements of SMC 
at CERN, and E143 and E155 at SLAC, but they are not 
statistically precise enough to detect a deviation of g2 
from its Wandzura-Wilczek part, as seen by the SLAC 
experiments. 
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Table 1. The spin-structure function xg2{x,Q^) and the virtual-photon asymmetry A2{x,Q^) of the proton in bins of {x,Q^), 
see text for details. Statistical and systematic uncertainties are presented separately 



bin 


(x) 


(0">,GeV^ 


X92 


± stat. 


± syst. 


A2 


± stat. 


± syst. 


1 


0.009 


0.38 


0.0799 





0521 





0182 


0.0257 


0.0163 


0.0057 


2 


0.018 


0.68 


0.0699 





0513 





0111 


0.0269 


0.0183 


0.0040 


3 


0.033 


0.89 


0.0450 





0326 





0215 


0.0278 


0.0165 


0.0109 


4 


0.039 


1.37 


-0.0047 





0652 





0080 


0.0033 


0.0275 


0.0035 


5 


0.044 


1.80 


0.3489 





1279 





0612 


0.1440 


0.0507 


0.0243 


6 


0.067 


1.09 


0.0044 





0421 





0097 


0.0190 


0.0346 


0.0085 


7 


0.069 


1.88 


0.0473 





0357 





0062 


0.0402 


0.0210 


0.0041 


8 


0.076 


2.79 


0.0202 





0674 





0323 


0.0225 


0.0342 


0.0164 


9 


0.116 


1.30 


-0.0094 





0506 





0081 


0.0266 


0.0603 


0.0111 


10 


0.118 


2.44 


0.0356 





0301 





0099 


0.0584 


0.0251 


0.0090 


11 


0.124 


4.04 


-0.0571 





0466 





0149 


-0.0137 


0.0311 


0.0102 


12 


0.182 


1.51 


-0.0758 





0642 





0230 


-0.0466 


0.1055 


0.0389 


13 


0.183 


3.01 


0.0121 





0324 





0038 


0.0707 


0.0375 


0.0074 


14 


0.187 


5.42 


-0.0334 





0440 





0041 


0.0143 


0.0392 


0.0052 


15 


0.282 


1.95 


0.0071 





0396 





0063 


0.1675 


0.0925 


0.0167 


16 


0.298 


3.99 


-0.0242 





0195 





0055 


0.0718 


0.0363 


0.0117 


17 


0.311 


7.58 


-0.0571 





0283 





0105 


0.0039 


0.0437 


0.0166 


18 


0.458 


2.83 


-0.0613 





0582 





0129 


0.0064 


0.2616 


0.0598 


19 


0.482 


4.31 


-0.0987 





0370 





0104 


-0.2064 


0.1704 


0.0500 


20 


0.484 


7.57 


-0.0362 





0183 





0045 


0.0421 


0.0744 


0.0206 


21 


0.630 


4.76 


0.2413 





1194 





0534 


3.0231 


1.3295 


0.5969 


22 


0.658 


6.79 


-0.0129 





0320 





0081 


0.1197 


0.4350 


0.1115 


23 


0.678 


10.35 


0.0076 





0160 





0025 


0.3672 


0.2551 


0.0419 



Table 2. The spin-structure function xg2 and the virtual-photon asymmetry A2 of the proton after evolving to common 
and averaging over in each x-hm (see text for details). Statistical and systematic uncertainties are presented separately 



bin 


X range 


{x) 


{Q">,GeV^ 


xg2 


zbstat 


zbsyst 


A2 


zbstat 


zbsyst 


1 


0.004 


- 0.014 


0.009 


0.38 


0.0794 


0.0520 


0.0153 


0.0256 


0.0162 


0.0049 


2 


0.014 


- 0.023 


0.018 


0.68 


0.0668 


0.0509 


0.0181 


0.0258 


0.0182 


0.0065 


3 


0.023 


- 0.050 


0.036 


1.08 


0.0456 


0.0262 


0.0157 


0.0261 


0.0121 


0.0074 


4 


0.050 


- 0.090 


0.069 


1.59 


0.0271 


0.0236 


0.0150 


0.0312 


0.0154 


0.0100 


5 


0.090 


- 0.150 


0.118 


2.07 


-0.0023 


0.0212 


0.0085 


0.0289 


0.0194 


0.0088 


6 


0.150 


- 0.220 


0.183 


2.51 


-0.0005 


0.0086 


0.0063 


0.0612 


0.0109 


0.0105 


7 


0.220 


- 0.400 


0.291 


3.23 


-0.0314 


0.0126 


0.0043 


0.0629 


0.0248 


0.0104 


8 


0.400 


- 0.600 


0.473 


4.62 


-0.0454 


0.0154 


0.0075 


0.0373 


0.0665 


0.0345 


9 


0.600 


- 0.900 


0.654 


7.06 


0.0107 


0.0177 


0.0073 


0.4275 


0.2316 


0.0970 



Table 3. Correlation matrix for xg2 in 9 a;-bins (as in Table [2]) 





1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


1.0000 


-0.1281 


-0.0038 


-0.0033 


-0.0017 


0.0005 


0.0000 


0.0000 


0.0000 


2 


-0.1281 


1.0000 


-0.1584 


-0.0083 


-0.0007 


0.0000 


0.0000 


0.0001 


0.0000 


3 


-0.0038 


-0.1584 


1.0000 


-0.1951 


-0.0281 


0.0077 


-0.0016 


0.0002 


0.0000 


4 


-0.0033 


-0.0083 


-0.1951 


1.0000 


-0.2885 


0.0312 


-0.0107 


0.0013 


-0.0005 


5 


-0.0017 


-0.0007 


-0.0281 


-0.2885 


1.0000 


-0.0102 


-0.0654 


0.0067 


-0.0018 


6 


0.0005 


0.0000 


0.0077 


0.0312 


-0.0102 


1.0000 


-0.1829 


0.0143 


-0.0055 


7 


0.0000 


0.0000 


-0.0016 


-0.0107 


-0.0654 


-0.1829 


1.0000 


-0.3539 


0.0926 


8 


0.0000 


0.0001 


0.0002 


0.0013 


0.0067 


0.0143 


-0.3539 


1.0000 


-0.3947 


9 


0.0000 


0.0000 


0.0000 


-0.0005 


-0.0018 


-0.0055 


0.0926 


-0.3947 


1.0000 



